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I.  INTRODUCTION 


For  many  years  it  has  been  observed  that  liquid-filled  projectiles  have  a 
proclivity  for  unusual  flight  behavior,  often  being  unstable  even  though  the 
same  projectile  with  a  solid  payload  is  stable.  The  basic  states  of  motion  of 
the  liquid  payload  in  a  cylinder  can  be  rationally  separated  into  solid-body 
rotation  and  spin-up,  the  transient  state  that  ends  with  the  fluid  in  solid- 
body  rotation.  The  determination,  according  to  linear  theory,  of  the  fre¬ 
quencies  and  decay  rates  of  the  waves  in  the  rotating  fluid  is  called  the 
"eigenvalue  problem."  It  was  solved  for  the  basic  state  of  solid-body  rota¬ 
tion  in  Reference  1.  The  purpose  of  this  report  is  to  present  the  theory  for 
the  basic  state  of  spin-up  and  some  representative  results. 

Some  attempts  to  solve  the  eigenvalue,  abbreviated  e.v.,  problem  have 
been  made  in  the  past.  Reddi2  attempted  to  solve  the  inviscid  perturbation 
equations  but  had  difficulty  with  a  singularity  in  these  equations,  related  to 
the  critical  layer  to  be  discussed  in  Chapter  IV.  A  critique  of  this  work  is 
given  by  Lynn.3  An  ad  hoc  method  to  solve  the  e.v.  problem  was  considered 
in  Reference  4.  The  actual  spin-up  velocity  profile  was  replaced  by  distri¬ 
butions  in  which  the  angular  velocity  was  constant  in  an  annular  layer  next  to 
the  cylinder  and  zero  elsewhere.  These  distributions  were  required  to  have 
either  the  same  angular  momentum  or  the  same  volume  of  rotating  fluid  as  the 
correct  spin-up  profile.  In  this  way  the  spin-up  problem  could  be  fitted  into 
the  formalism  of  the  Stewartson  theory.5  This  ad  hoc  approach  was  not 
successful,  and  the  e.v.'s  disagreed  with  results  from  the  present  theory.  A 
rational  approach  to  the  e.v.  problem  was  given  by  Lynn3  for  a  basic  flow 
designated  by  Vw  in  Section  II.  The  effect  of  the  critical  layer  singularity 

was  discussed  and  the  real  part  of  the  e.v.  calculated  for  axisymmetric 
inviscid  disturbances,  for  which  no  critical  layer  exists. 


1.  C.  U.  Kitchens,  Jr. ,  V.  Gerber,  and  R.  Sedney,  "Oscillations  of  a  Liquid 
in  a  Rotating  Cylinder:  Paid:  I.  Solid  Body  Rotation ARBRL-TR-02081 , 
June  19? 8  (At)  A057Y59)  , 

3.  M.  M,  Reddi ,  " On  the  Rigenvalueo  of  Couatte  Flow  in  a  Fully-Filled 
Cylindrical  Container Franklin  Institute  Reseai'eh  Laboratory, 
Philadelphia ,  PA,  Report  So.  F-B2394,  1987. 

3.  T.  M.  Lynn ,  "Free  Oscillations  of  a  Liquid  During  Spin-Up,"  BRL  Report  So. 
1663,  August  1973  (AD  769710). 

4.  Fngineenng  Design  Handbook,  Liquid-Filled  Projectile  Design,  A>4C  Pamphlet 
So.  706.165,  United  States  Amy  Hateriel  Comtand,  Washington,  D.C.,  April 
1969  (AD  853719). 

5.  K.  Staeartson,  '* On  the  Stability  of  a  Spinning  Top  Containing  Liquid," 
Journal  of  Fluid  Mechanics ,  Vol.  S,  Part  4,  1959. 
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I.i  this  report  the  linear  e.v.  problem  is  solved  for  viscous  pertur¬ 
bations,  axi symmetric  or  not,  on  a  basic  flow  which  is  the  spin-up  core  flow 
including  diffusion.  The  modal  analysis  used  allows  slip  on  the  cylinder 
endwalls.  The  perturbation  equations  have  a  singularity  at  the  axis  of  the 
cylinder  that  is  handled  analytically.  Because  the  Reynolds  number  is  large, 
an  integration  scheme  with  orthonormalization  is  used  that  maintains  linearly 
independent  solutions.  The  determination  of  the  eigenfunctions  through  the 
critical  layer  is  accomplished  numerically. 

The  significance  of  spin-up  effects  on  a  projectile  flight  can  be  esti¬ 
mated  in  an  order  of  magnitude  sense.  Let  a  and  c  be  the  radius  and  half¬ 
height  of  the  cylinder,  ft  the  spin  in  rad/sec  of  the  projectile  (say  at  the 
muzzle),  v  the  kinematic  viscosity  of  the  liquid,  and 

2  2 
Re  =  a  a  /v,  E  -  \>/'ft  c  , 

the  Reynolds  number  and  Ekman  number.  Note  that  E  is  sometimes  defined  using 
the  height,  2c,  rather  than  c,  and  that  for  historical  reasons  we  shall  use  Re 
in  this  report  rather  than  the  more  conventional  E.  If  the  Ekman  layers, 
i.e.,  endwall  boundary  layers,  are  laminar,  the  characteristic  time  for  spin-up 
is 

ts  «  (2c/a)Re1/2/n 
=  2/E1/2ft  (sec). 

We  also  use  a  nondimensional  spin-up  time  t$  a  ft  For  Re  >  105, 
approximately,  the  Ekman  layer  may  be  turbulent,  in  which  case  the  charac¬ 
teristic  spin-up  time  can  be  estimated  by  £^=(28.6  c/a)  Re^/ft. 

Basically,  these  are  time  scales  for  the  spin-up  process  and  do  not  neces¬ 
sarily  give  a  measure  of  how  close  the  spin-up  process  is  to  solid-body 
■citation.  This  can  be  measured  in  several  ways,  e.g.,  time  for  the  angular 
momentum,  e.v.,  or  other  property  to  attain  a  value  within  10%,  say,  of  the 
solid  body  value.  These  measures  require  detailed  calculations  and  so  are  not 
very  convenient.  A  rule  of  thumb  is  that  solid-body  rotation  is  reached  at 

about  4  t  after  an  impulsive  start  of  the  cylinder. 

The  characteristic  spin-up  times  can  be  compared  with  tfl  ,  the  time  of 

flight  of  the  projectile.  If  t$  or  t  «  tfl,  spin-up  effects  can  be 

disregarded  and  solid-body  rotation  can  be  assumed.  If  tc  or  t$t  a  t^/10, 

or  larger,  spin-up  effects  probably  have  to  be  considered.  To  put  these 
estimates  in  perspective,  consider  the  parameters  for  two  of  several  recent 
firings  of  155aw  projectiles  conducted  by  Or  W.  P.  D'Amico,  BRL.  For  one  of 
these  c/a  -  3.120,  Re  *  4  x  104,  ft  784  rad/sec,  giving  t  e  1.65  sec;  for 

another  case  c/a  =  5.200,  Re  =  2  x  10*>,  ft  =  754  rad/sec  giving,  £  =  3.58 

sec.  These  are  referred  to  as  Cases  1  and  2,  respectively.  For  both  of  these 

tfj  *=  40  sec.  and,  from  the  e.v.  histories,  spin-up  effects  are  significant. 
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A  question  of  more  practical  concern  is:  To  what  extent  does  the  spin-up 
process  affect  the  moment  exerted  by  the  liquid  on  the  projectile  and  its  yaw 
growth  rate?  The  answer  requires  solution  of  the  forced  oscillation  problem 
for  the  fluid.  The  e.v.  problem  is  the  free  oscillation  problem,  and  its 
role,  vis-a'-vis  the  forced  oscillation  problem,  is  the  same  here  as  in  any 
vibrating  system.  For  example,  the  e.v.  and  eigenfunctions  can  be  used  to 
expand  the  solution  of  the  forced  oscillation  problem.  The  e.v.  can  also  be 
used  to  estimate  when  the  projectile  might  be  unstable.  For  inviscid 
perturbations  on  solid-body  rotation,  Stewartson5  showed  that  if  the  projec¬ 
tile  nutational  frequency  is  within  a  certain  band  of  the  eigenfrequencies, 
the  moment  has  a  resonant  response  and  the  projectile  will  be  unstable.  If 
this  idea  is  extended  to  viscous  perturbations  on  spin-up,  the  real  part  of 
the  e.v.  would  be  compared  with  the  nutational  frequency.  The  comparison  is, 
at  best,  a  necessary  condition  for  projectile  instability.  Confidence  in 
application  of  this  theory  must  be  obtained  by  other  means.  This  has  been 
provided  by  experiments11  and  numerical  simulations.7 


II.  SPIN-UP  BASIC  FLOW 

Before  the  spin-up  e.v.  problem  is  considered, the  basic  flow  requires 
some  discussion.  Only  those  aspects  of  it  germane  to  the  e.v.  problem  will  be 
discussed.  Our  interest  is  in  spin-up  from  rest  which  is  inherently  a  non¬ 
linear  problem.  Spin-up  from  one  state  of  solid-body  rotation  to  another  can 
be  treated  as  a  linear  problem  if  the  change  in  rotation  rate  is  small.  Both 
types  are  discussed  by  Greenspan**  and  Benton  and  Clark.' 3 

Consider  the  motion  of  a  fluid  which  fills  a  cylinder,  initially  at  rest, 
which  is  "rapidly"  brought  to  a  constant  angular  velocity.  The  fluid  motion 
is  axi symmetric  and  time  dependent.  There  are  many  variations  on  this  problem 
which  are  of  practical  interest,  e.g.,  the  fluid  only  partially  filling  the 
cylinder,  but  they  will  not  lie  considered  here.  For  our  application  the 
angular  acceleration  of  the  cylinder  is  large  during  the  time  when  the  pro¬ 
jectile  accelerates  in  the  gun  tube;  the  small  deceleration  in  flight  is 
neglected.  For  artillery  projectiles  this  acceleration  time  is  typically 

0.020  sec.  which  is  small  compared  to  ;  comparison  of  angular  acceleration 


6.  S.  Stergiopouloa ,  "An  Experimental  Study  of  Inertial  Raves  in  a  Fluid 
Contained  in  a  Rotating  Cylindrical  Cavity  During  Spin-Up  Frm  Root/'  Ph.D. 
Thesis ,  York  University ,  Toronto ,  Ontario,  February  1982. 

7.  R.  Sedney,  N.  Gerber,  art!  J.H.  dartos,  "Oscillations  of  a  Liquid  in  a 
Rotating  Cylinder,"  4/44  Preprint  82-0296,  20th  Aero&[Xice  Sciences  Heating, 
Orlando,  FL,  January  1982. 

8.  H.  P.  Greensfxw,  The  'Ah cony  of  Rotating  Fluids,  Cambridge  University  Press , 
London  and  N&e  York,  1988. 

9.  E.  R.  Renton  and  A.  Clark,  Jr.,  "Spin-lip,"  article  in  Annual  Re  id  on  of 
Fluid  Hschanics,  Vol.  8,  4«nual  Reoieus,  Inc.,  Palo  Alto,CA,  id74. 
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time  with  spin-up  time  is  also  required.  The  assumption  of  an  Impulsive  start 
for  the  cylinder  is  a  good  approximation  for  the  projectile  problem;  it  may 
not  be  in  a  laboratory  apparatus.  For  an  impulsive  start  our  model  requires 
separate  treatment  of  the  initial  conditions. 

The  notation  here  will  be  the  same  as  in  Reference  1.  Lengths,  veloci- 

2  2 

ties,  pressure,  and  time  are  made  nondimensional  by  a,  an,  pa  a  ,  and 

a  *, respectively,  where  p  is  the  liquid  density  and  n  is  the  constant  spin 
rate  of  the  cylinder.  In  the  Inertial  frame  cylindrical  coordinates  r,  e,  2 
are  used,  with  the  origin  of  z  at  the  center  of  the  cylinder,  and  velocities 
U,  V,  W,  respectively.  Dimensionless  time  Is  t.  Note  that  in  Reference  1  the 
origin  of  z  is  at  the  bottom  endwall.  Derivatives  are  indicated  by  sub¬ 
scripts. 

The  basic  work  on  spin-up  from  rest  was  done  by  Wedemeyer.10  In  this 
remarkable  paper  he  shows  that  the  flow  must  be  divided  into  two  regions:  the 
Ekman  layers  at  the  endwalls,  which  can  be  treated  in  a  quasi-steady  fashion, 
and  the  rest  of  the  flow,  called  the  core  flow.  Wedemeyer  did  not  point  out 
that  an  additional  boundary  layer,  a  Stewartson  layer,  is  squired  at  the 
cylinder  wall.  Starting  with  the  Navier-Stokes  equations  for  axisymmetric 
flow: 


u*  +  u*  u* 

t  r 

+  w*  u* 

2 

? 

-  V*£/r  = 

-P*  + 
r 

Re'1 

(V2U*  -  U*/r2) 

(2.1a) 

V*  ♦  U*  V*  ♦ 
t  r 

W*  V*  ♦ 

U*  V*/r  ■ 

Re'1 

<V 

-  V*/r2) 

(2.1b) 

W* 

t 

♦  u*  w* 

r 

+  w*  u*  ■ 

l 

-PJ  ♦ 

Re'1 

v2  w* 

(2.1c) 

<rU*) 

r*  r“I‘ 

0 

(2. Id) 

v2(!  ■ 

<>rr 

♦  U/r)  ()p*  ()„. 

U.le) 

where  the  asterisk  indicates  the  exact  solution,  he  used  physical  order-of- 
magnitude  arguments  to  simplify  them  in  the  core  flow.  The  approximations 
reduce  the  three  momentum  equations,  (2.1,  a,  b,  c),  to 

Vt  ♦  U  (Vp  ♦  V/r)  ■  Re'1  [Vrr  *  <V/r)p3  (2.2) 


10,  S,  H.  Wedewyer,  "The  Unsteady  Flou  Within  a  Spinning  Cylinder,"  USA  BRL 
Report  No,  1825 ,  October  1963  (AD  431346).  Also,  Journal  of  Fluid 
Mechanics ,  Vol*  20,  Part  3,  1964,  pp.  383-399. 
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and 


Uz  -  V2  =  P2  =  0,  (2.3) 

where  the  superscript  is  dropped  for  this  approximate  solution.  For  Re  +  ® 
Wedemeyer  proposed  neglecting  the  diffusion  terms  in  (2.2)  and  arrived  at 

v  +  Uw  (V  ♦  Vw/r)  =  0,  (2.4) 

where  the  sub  w  indicates  this  approximation.  Wedemeyer  used  (2.4)  rather 
than  (2.2)  when  he  applied  his  model. 

Although  Wedemeyer  showed  the  crucial  importance  of  the  Ekman  layers  to 
the  spin-up  process,  his  model  did  not  require  a  solution  for  the  flow  in 
these  layers.  He  did  discuss  certain  properties  of  the  solution  which  were 
required  in  his  model  for  the  core  flow;  exclusion  of  this  solution  has 
important  consequences  for  the  e.v.  problem. 

A  more  formal  approach  to  this  problem  was  given  by  Greenspan.8  In  his 
treatment,  lengths  are  made  dimensionless  by  2c  rather  than  a,  which  is  more 

natural  for  the  spin-up  process;  time  is  made  dimensionless  by  £  so  that  the 

1 12  S 

new  time  is  t*  *  t/t  =  t/(2c/a)  ReA/  ;  an  expansion  in  the  small  parameter 

*5 

<l/ts)  -  (a/2c)  Re‘i/2  (2.5) 

is  assumed.  The  form  of  this  expansion 
U*  «  (l/ts)  Uj  ♦  .  .  . 

V*  -  VQ  ♦  (i/t$)  ♦  .  .  .  (2.6) 

W*  •  (l/ts)  Uj  «•  .  .  . 

p*  3  p0  *  U/ts)  *  .  .  . 

follows  from  the  facts  that  in  the  core  flow  U*  and  W*  are  0(  1/t s )  and  V*  is 

0(1),  if  1/t s  is  small.  See  Reference  8  for  details.  The  independent 

variables  are  r,  2,  t'.  Substituting  the  expansion  into  (2.1)  after  trans¬ 
forming  to  t'  yields 
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(2.7) 


V  +  U.  V  +  V  /r  =  0 
0.  ,  1  '  0  o  ' 

t  r 


V02/r  « 


r 


These  are  the  same  as  (2.3)  and  (2.4)  since  V  =  V  ,  U  =  (1/t  )  U.  ,  and 

t1  =  t/ts,  Although  the  formalism  of  matched  asymptotic  expansions  was  not 

used,  the  first  two  terms  of  the  outer  expansion  are  apparently  given  by 
(2.6).  The  Ekman  layer  solution  would  be  obtained  from  the  inner  expansion 
which  was  not  considered  in  Reference  8. 

To  solve  (2.2),  (2.4),  or  (2.7)  a  relationship  between  U  and  V  is 
necessary.  Wedemeyer  used  the  facts  that  the  Ekman  layers  are  steady  after 
one  revolution  and  that  the  radial  mass  flux  in  the  core  flow  must  be  balanced 
by  that  in  the  Ekman  layers  to  obtain  some  conditions  on  the  U,  V  relation¬ 
ship.  At  this  point  he  was  forced  to  take  a  phenomenological  approach.  The 
relationship  is  known  at  t  ♦  0  and  t  *  and  he  proposed  a  linear  inter¬ 
polation  between  them  to  obtain  an  approximate  >*elationship  for  any  t.  He 
tested  this  idea  in  sane  other  problems  where  the  solution  was  known  ami 
decided  it  was  satisfactory.  Some  confusion  has  appeared  in  the  literature 
because  this  step  was  misinterpreted;  this  is  discussed  in  a  separate 
report.  The  result  is 


U  ■  (V  -  r) 


k t  *  n  (a/c)  Re~^ 


(2.8) 


£  1/2 


for  laminar  Ekman  layers.  Wedemeyer  proposed  *  =  0.443  but  Greenspan  sug¬ 
gested  c  *  0.5;  the  latter  often  gives  results  in  better  agreement  with 
numerical  solutions  to  the  Navier-Stokes  equations.  Otherrelationships  have 
been  proposed,  but  they  will  not  be  discussed  here.  For  turbulent  Ekman 
1 ayers 
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(2.9) 


&,v-.-’JWWpi^  r.'V-P^^K.'T^vp-j 


U  =  -kt  (r-V)8/5 
kt=  0.035  (a/c)  Re_1/5 

’  ^st* 


where  t  ^  is  the  nondimensional  ,  turbulent  spin-up  time..  The  core  flow  is 

assumed  to  be  laminar  so  that  turbulent  stresses  are  not  introduced  in  the 
right-hand  side  of  (2.2), 

Using  (2.8),  (2.4)  can  be  solved  explicitly: 


2k  t  2k  t  -k  t 

V  =  (re  -  l/r)/(e  -  1)  for  r  >  e 

W 


~kil 

for  r  <  e  , 


(2.10) 


~kl- 

so  that  r  =  e  separates  rotating  and  nonrotating  fluid.  Although  V 

W 

is  continuous  there,  a  discontinuity  in  shear,  i.e.,  V  ,  exists  there.  Using 

wr 

(2.9)  a  numerical  integration  is  necessary  to  obtain  V^;  but  the  character  of 

the  solution  is  the  same  as  when  (2.8)  is  used.  The  radial  velocity  is 
obtained  from  (2.8)  or  (2.9)  and  W  is  obtained  from  the  continuity  equation 


W  =  -(z/r)  (rU)p 

At  r  -  1,  W  #0;  thus,  the  Stewartson  layer  should  be  included  at  r  =  1,  but 
this  has  yet  to  be  done.  Also  W  *  0  at  the  endwalls  z  ®  i  c/a. 

Wedemeyer  pointed  out  that  the  corner  in  the  solution  to  (2.4)  would  be 
smoothed  out  if  the  diifusion  terms  were  retained,  as  in  (2.2).  But  this  is  a 
nonlinear  second  order  equation  and  must  be  integrated  by  finite  difference 
methods.  This  can  be  done  with  some  standard  techniques  for  diffusion  equa¬ 
tions,  using  (2.Q)  or  (2.9).  Special  treatment  is  needed  near  the  point 
r  -  1,  t  =  0  because,  for  an  impulsive  start,  a  discontinuity  in  the  boundary 
conditions  exists  at  that  point.  In  our  work  a  local,  analytic  solution  was 
derived  to  resolve  the  discontinuity.  In  most  of  our  e.v.  calculations  we 
used  the  V  from  the  numerical  solution  of  (2,2)  with  either  (2.8)  or  (2.9), 
but  other  options  are  also  used.  For  laminar  Ekman  layers,  the  spin-up 
velocity  profile 

V  3  f  (r,  kjt,  kt  Re); 


13 


■;fe  Uft 


for  the  turbulent  case  is  replaced  by  k^„  These  functional  forms  follow 

directly  from  (2.2)  and  (2.8)  or  (2.9)  and  represent  the  most  efficient  way  to 
display  the  results. 

Some  examples  of  the  solution  of  (2.2)  will  De  giver,  simply  to  illustrate 
the  form  of  V  which  must  be  perturbed  in  the  e.v.  problem.  The  parameters  of 
the  two  cases  of  155mm  projectiles,  given  in  the  introduction,  will  be  used. 
Case  1,  with  Re  =  4  x  104,  should  have  laminar  Ekman  layers.  The  solution  to 
(2.2)  with  (2.8)  giyes  the  V  vs  r  shown  in  Figure  1  for  3  values  of  t.  One  of 
them  is  the  nondimensional  spin-up  time  t  =  1245.  For  t  =  4800,  approxi¬ 
mately  4t  ,  solid-body  rotation  is  achieved,  essentially.  Case  2,  with  Re 
=  2  x  IQ6,  should  have  turbulent  Ekman  layers.  The  solution  to  (2.2)  with 

(2.9)  gives  the  V  vs  -  shown  in  Figure  2  for  3  values  of  t  including 

t$t  =  2700.  For  c  =  10,000,  approximately  4  t  ,  solid-body  rotation  has  not 

been  achieved  to  the  same  degree  as  in  Figure  1.  The  profile  for  t  =  1300  is 
relatively  steep;  early  time  profiles  such  as  this  require  more  care  in  the 
e.v.  calculation. 

Some  results  from  the  Wedemeyer  spin-up  model  have  been  compared  with 
experiments;  these  tend  to  validate  toe  model.  Another  way  to  validate  the 
model  is  to  compare  its  results  with  solutions  to  the  Navier-Stokes  equa¬ 
tions.  The  firJte  difference  solution  presented  by  Kitchens11  was  used  for 
this  purpose. 

The  Vs  from  both  solutions  were  compared  over  a  significant  range  of 
parameters;  in  general,  the  differences  were  small  for  k  =  0.5  and  slightly 
larger  for  <  =  0.443.  Comparisons  are  given  in  a  separate  report.  Some 
limited  comparisons  are  shown  in  the  following  discussion  which  actually  has  a 
different  purpose,  viz.,  to  show  the  flow  in  the  Ekman  layer  and  its  transi¬ 
tion  to  the  core  flow. 

To  the  authors'  knowledge  no  solution  for  the  flow  in  the  Ekman  layers 
during  spin-up  has  been  given  other  than  finite  difference  solutions.  In 
References  8  and  10  certain  properties  of  these  were  used  to  determine  the 
core  flow.  The  finite  difference  solutions  to  the  Navier-Stokes  equations  by 
Kitchens11  provide  the  complete  flow.  In  Figure  3  results  from  this  method 
are  shown  for  c/a  =  1,  Re  *  104,  t  =  200  at  t  °  62.5  and  250.  For  three 

values  of  r,  V*  vs  z  +  c/a  is  shown.  According  to  Wedemeyer's  spin-up  model  v 
is  independent  of  z;  this  is  verified  by  the  calculation  for  z  +  c/a  .05  at 
t  =  62.5  and  for  z  +  c/a  >  .025  at  t  =  250.  These  values  of  z  +  c/a  can  be 
taken  as  one  measure  of  Ekman  layer  thickness.  The  values  of  V  from  the 
Wedemeyer  model  are  also  indicated;  the  differences  from  the  results  of  the 
Navier-Stokes  equation  are  relatively  large  at  t  =  62.5  and  essentially  zero 
at  t  =  250,  Significant  undershoots  exist  in  the  V  pi'ofile  at  t  =>  62.5  but 


n.  C ,  </.  Kitchens,  Jr.,  "Navier-Stokes  Solutions  for  Spin-Up  in  a  Filled 
Cylinder,"  A IAA  Journal .  Vol.  18,  No.  8,  August  1980,  pp .  929-934 .  Also 
Teohnioal  Report  ARBRL-TR-Q2193 ,  September  1979  (AD  A07711S). 
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only  mild  ones  at  t  =  250.  An  example  of  U  vs  z  +  c/a  is  given  in  Reference 
11.  The  azimuthal  component  of  the  basic  flow  which  must  be  perturbed  to 
obtain  the  e.v.  is  represented  in  Figure  3. 

If  only  the  core  flow  were  perturbed,  it  would  not  be  possible  to  satisfy 
the  boundary  conditions  on  the  endwalls.  Figure  3  also  shows  that  a  large 
error  in  azimuthal  velocity  is  made  in  the  Ekman  layers  if  the  core  flow  is 
extended  to  the  endwalls.  This  is,  in  fact,  what  is  done  in  the  e.v,  problem 
formulated  in  this  report.  Nevertheless,  the  wave  frequency  calculated  in  this 
way  agrees  fairly  well  with  some  experimental  results;  however,  the  calculated 
wave  damping  has  a  substantial  error.  This  is  expected  since  the  dissipation 
in  the  Ekman  layers  should  affect  the  wave  damping. 


III.  THE  PERTURBED  FLOW 

The  procedure  for  obtaining  the  equations  governing  the  perturbed  flow  is 
a  standard  one.  It  starts  from  the  Navier-Stokes  equations  for  3-D  flow: 

U*  U  U  +  (V/r)  U,  +  W  U  -  V2/r  =  -P,.  +  Re_1(v2  U  -  U/r2  -  2V Jr?) 
v  r  6  z  r  $ 

V.  +  UV  +  (V/r)  Vft  +  W  V,  U  V/r  -  -Pfl/r  +  Re”1  (v2  J  -  V/r2  +  2U  Jr2) 
t  r  u  z  □ 

W  +  II  W  +  (V/r)  W  +  W  W  =  -P  +  Re"V  W  (3.1) 

t  r  e  z  z  v  *  ' 

(r  U)r  +  VQ  +  r  W2  =  0 

)  ■  V,  (  )  +  ( 1/r2)  (  )  oe 


The  velocity  components  and  pressure  are  then  expressed  as  the  sum  of  the 
basic  and  the  perturbation: 


U  »  U*  (r,z,t)  +  u*  (r,e,z,t)  (3.2) 


and  similarly  for  V,  W,  and  P.  For  our  application  nonaxi symmetric 
disturbances  must  be  considered,  hence  the  dependence  of  u1  on  0. 

The  details  of  the  derivation  will  be  given  only  for  the  radial  momentum 
equation,  the  first  of  (3.1),  The  representation  of  the  velocity  components 
and  pressure  as  a  basic  flow  plus  a  perturbation,  i.e.,  (3.2)  and  the  similar 

expressions  for  V,  W,  and  P,  are  substituted  into  (3,1).  The  terms  are 
separated  into  zeroth  order  terms  which  are  independent  of  the  u* ,  v* ,  w1 ,  p' , 
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first  order  terms  which  are  linear  in  u* ,  v1,  w‘ ,  p1 ,  etc.  The  zeroth  order 
terms  give  exactly  (2.1),  by  construction.  The  first  order  terms  for  the 
radial  momentum  equation,  the  first  of  (3.1),  give 


u‘t  +  U*u‘r  +  u'  U*r  +  (V*/r)  u 1 q  +  W*  u‘z+  w1  U*2  -  (2V*v')/r  =  (3.3) 


-p'  +  Re"1  [T2  u'  -  u'/r2  -  2v'  /r2]. 

One  term  does  not  appear  because  U*  is  independent  of  e. 

Uo  to  this  point  the  basic  flow  is  an  exact  solution  of  the  Navier-Stokes 
equations,  (2,1),  which  would  have  to  be  a  finite  difference  solution.  It  is 
impractical  to  use  such  a  solution  for  the  coefficients  in  (3.3).  Therefore, 
an  approximation  to  it,  the  Wedemeyer  model  for  the  core  flow,  is  used.  This 
leaves  open  the  matter  of  treating  the  Ekman  layer  flow  and  the  attendant 
difficulties  discussed  ct  the  end  of  the  last  section.  A  correction  is 
"equirad  to  coui.u  for  the  Ekman  layers;  this  will  not  be  done  here. 

The  next  step  is  to  estimate  the  order  of  magnitude  of  the  terms  in 
(3.3).  They  are  not  al i  the  same  order  of  magnitude,  and  this  step  leads  to  a 
significant  s .mplification.  The  expansion  (2.6)  is  convenient  for  this  pur¬ 
pose  although  the  approach  t„xen  in  deriving  (2.2)  was  actually  used.  It  is 
assumed  that  al1  pertui bations ,  the  primed  quantities,  and  their  derivatives 
are  0(1).  in  particular  u'  -  0(1)  which  means  that  the  dimensional  ti>«e 

scale  for  the  perturbations  s  Ohi  )  and  is  therefore  small  compared  to  t  ; 

the  nondimensional  time  scale  for  the  perturbations  is  0(1)  which  is  small 
compared  to  This  estimation  is  used  below  to  make  another  important. 

simplification.  It  follows  fro, a  (2.6),  or  the  equivalent  using  Wedemeyer's 
approach,  that  all  coefficients  in  (3,3)  '•oniaining  U*  and  W*  and  their 
derivatives  with  respect  to  r  and  z  are  0(l/t&);  all  those  containing  V*  and 

V*r  are  0(1).  Neglecting  the  0(l/ts)  terms  on  the  left-hand  side  of  (3.3) 

deletes  all  terms  that  contain  U*  and  W*.  Only  the  term.,  containing  V*  and 
those  not  involving  the  basic  flow  remain.  For  V*  we  must  substitute  V  of 
(2.2).  Applying  this  same  process  to  remaining  equations  of  (3.1)  yields 
the  following  set  of  perturbation  equations; 

u't  +  (V/r)  u'g  -  2Vv'/i  *  -p‘r  +  R''~1  (v2u*  -  u‘/r2  -  2v*g/r2)  (3.*a) 

v‘t  +  [Vr  ♦  (V/r)]  u‘  ♦  (V/r)  v’  Q  -  -p’  Jr  *  Re  1  (7  V  (3.4b) 

-  v‘/r2  +  2u‘  fr?) 
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(3.4c) 


w'  +  (V/r)  w*  =  -p*  +  Re"1  v2  w1 

L  o  Z 

(ru')r  +  vjj  +  r  =  0.  (3.4d) 

These  are  the  equations  that  govern  viscous  perturbations  of  the  core  flow  but 
not  the  Ekman  layers. 

The  presence  of  the  viscous  terms  on  the  right-hand  side  of  (3.4a,  b,  c) 

requires  discussion.  It  is  helpful  to  refer  back  to  the  spin-up  basic  flow 

model.  A  formal,  rational  expansion  such  as  (2.6)  yields  (2.7)  in  which  no 

viscous  diffusion  terms  appear,  but  Wedemeyer's  approach  did  have  these  terms, 

as  in  (2.2).  Without  them  the  predictions  of  the  model  would  not  agree  nearly 

as  well  with  experimental  or  other  numerical  results  for  a  range  of  r  and  t 

(See  Reference  12  for  a  discussion  of  this.)  As  mentioned  before,  at  the 

discontinuity  in  V  the  diffusion  terms  are  required.  They  could  be 
w 

r 

included  in  a  "local"  sense  as  Venezian13  did  or  in  a  "global"  sense  as 
Uedemeyer  did.  The  situation  with  (3.4)  is  analogous.  A  formal,  rational 
expansion  would  yield  perturbation  equations  of  the  same  form  as  (3.4)  but 

with  Re-1  =  0,  i.e.,  the  inviscid  equations.  These  are  adequate  except  for 
two  regions  where  viscous  effects  are  important.  These  are  the  boundary  layer 
at  r  -  1  and  the  critical  layer,  to  be  discussed  below.  Here  also  the  viscous 
terms  could  be  included  in  a  local  sense  in  both  of  these  regions  or  in  a 

global  sense  as  we  have  done  by  retaining  the  Re"1  terms  in  (3.4). 

There  is  another  restriction  on  these  equations  arising  from  the  fact 
that  Wedemeyer's  model  for  the  core  flow  is  not  correct  at  r  =  1,  as  discussed 
in  Section  II.  Specifically, W  *  0  at  r  =  1.  One  can  either  say  that  the 
model  is  valid  only  outside  a  Stewartson  layer  at  r  =  1  or  such  a  layer  could 
be  incorporated.  In  the  second  alternative  an  additional  term  would  appear  in 

(3.4c)  arising  from  the  term  0  in  the  z-momentum  equation  of  (3.1).  The 

first  alternative  will  be  adopted  here,  but  the  presence  of  the  Stewartson 
layer  is  neglected  just  as  the  presence  of  the  Ekman  layer  is.  The  fact 
that  W  *  0  at  z  e  t  c/a  is  part  of  the  overall  problem  of  neglecting  the  Ekman 
layer  in  the  basic  flow. 

The  final  step  is  to  consider  the  applicability  of  the  quasi-steady 
assumption  to  (3.4).  If  V(r,tl  varies  significantly  with  t  over  the  time 
scale  0(1)  appropriate  to  (3.4),  then  several  coefficients  depend  on  both  r 


12.  W,  8.  Watkins  and  R.  G .  Hussey,  "Spin-Up  From  Rest:  Limitations  of  the 

Wedemeyer  Model,"  The  Physios  of  Fluids,  Vol.  16,  No.  9,  September  1973, 
pp.  1530-1531.  ~  ~ 

13.  G.  Venenian,  Topics  in  Paean  Engineering ,  Vol .  2,  pp.  87-96,  Gulf 
Publishing  Co^,  Houston,  Texas,  lS7o. 
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and  t  and  separation  of  variables  could  not  be  used  to  solve  (3,4).  As 
Dointed  out  above  the  time  scale  for  V,  the  spin-up  time,  is  0(t$ )  which  is 

large  compared  to  0(1).  Therefore,  over  the  time  scale  for  the  perturbations 

V  does  not  change  appreciably.  Thus,  t  can  be  regarded  as  a  parameter  in  the 
solution  of  (3.4);  i.e.,  we  can  use  the  quasi-steady  approximation.  A 
possible  exception  to  this  might  exist  for  t  0  and  an  impulsive  start. 

In  Lynn's  work3  the  same  equations  as  (3.4)  appeared  except  that  his 
equations  have  V  in  place  of  V.  The  basic  flow  used  by  Lynn  was  (0,  Vw,  0, 

P  )  which  is  an  exact  solution  to  the  steady  Navier-Stokes  equations;  this 

basic  flow  would  not  apply  in  the  two  wall  layers.  The  philosophy  behind  his 
approach  is  different  from  ours  even  though  the  equations  are  the  same  if 

V  =  V  .  Lynn  solved  only  the  inviscid  equations.  Re"*  =  0,  for  axisymmetric 

n 

perturbations. 

From  the  discussion  in  the  introduction,  it  is  clear  that  wave-type 
solutions  to  (3.4)  are  required.  Although  there  may  be  several  approaches  to 
obtaining  these,  only  one  besides  that  actually  employed  will  be  mentioned. 

If  the  sinusoidal  variation  in  8  and  t  of  (u* ,v‘ ,w* ,p' )  Is  factored  out,  a  set 
of  linear  partial  differential  equations  In  r  and  2  is  obtained.  These  might 
have  to  be  solved  numerically.  This  approach  would  be  especially  useful  if 
the  Ekman  layer  effects  were  Included.  No  serious  attempt  has  been  made  to 
use  this  approach.  Here  a  separation  of  variables  is  employed  to  obtain  modal 
type  solutions. 


IV.  THE  MODAL  SOLUTION  AND  EIGENVALUE  PROBELM 

In  this  section  it  is  convenient  to  change  the  origin  of  the  z-coordlnate 
to  the  bottom  endwall  which  was  the  convention  in  Reference  1.  Let  that 
coordinate  be  z'  so  that  z‘  a  z  +  (c/a);  z*  will  be  used  Instead  of  z. 

It  is  assumed  that  the  perturbation  can  be  expressed  as  a  superposition 
of  modes  or  a  triple  Fourier  expansion  in  0,  z',  and  t  with  coefficients 
functions  of  r.  It  is  convenient  to  use  complex  notation  and  express  the 
perturbation  as 


u'  «  Real  |u  (r)  cos  K  z'  exp  [1  (Ct  -  ms)]  |  (4.1a) 


v'  =»  Real  |v  (r)  cos  Kz'  exp  [i  (Ct  -  mo)]  j  (4.1b) 


w‘  =  Real  jw  (r)  sin  Kz'  exp  [i  (Ct  -  mo)]  j  \4.lc) 
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(4. Id) 


p'  =  Real  jp  (r)  cos  Kz‘  exp  [i  (Ct  -  me)]  j 
K  =  kna/2c. 

The  complex  quantities  u,  v,  w,  and  p  are  solutions  of  the  ordinary  differ¬ 
ential  equations  obtained  by  substituting  these  forms  into  (3.4).  The 
integers  m  =  0,  +1,...  and  k  =  1,  2,...  are  the  azimuthal  and  axial  wave 
numbers,  respectively.  For  axisymmetric  perturbations  m  =  0;  m  =  1  is  the 
value  relevant  to  the  projectile  problem,  but  we  shall  not  specify  m  at  this 
point.  The  nondimensional  complex  constant  C  =  +  i  Cj  is  the  eigenvalue 

of  the  system.  The  dimensional  wave  frequency  is  ft  and  the  decay  rate  or 

damping  is  Cj  ft. 

For  the  free  oscillation  problem  the  boundary  conditions  at  the  endwalls, 
z1  =  0,  2c/a,  are  u*  =  v'  =  w1  =  0  if  the  complete  flow  is  being  perturbed. 
Since  we  are  perturbing  only  the  core  flow, the  boundary  conditions  are  not  the. 
same.  From  (4.1)  w'  =  0  at  the  endwalls  but  not  u'  and  v1.  Evidently, the 
form  (4.1)  represents  the  first  term  of  an  outer  expansion  of  the 
perturbation.  The  inner  expansion  must  be  introduced  to  satisfy  the  no-slip 
boundary  conditions. but  for  this  the  inner  expansion  of  the  basic  flow,  the 
Ekman  layers,  must  be  perturbed.  In  lieu  of  this,  the  form  (4.1)  is  used 
without  endwall  correction;  i.e,,  there  is  slip  at  the  endwalls. 

The  no-slip  boundary  conditions  at  r  =  1  are  satisfied  by  requiring 


u  =  v  =  w  a  0 


at  r  =  1. 


(4.2) 


Three  more  boundary  conditions  are  necessary.  For  a  filled  cylinder  the 
boundary  conditions  at  r  =  0  depend  on  m;  they  are  presented  in  Reference  1 
and  are  simply  restated  here: 


ill  =  0 

u  =  v  =  dw/dr  =  0 

m  a  1 

u-iv  =  w  =  p  a  0 

(4.3) 

m  >  1 

U  a  y  =  w  ■  0 

Applying  the  two  differential  operators  common  to  (3.4a,b,c)  to  the 
generic  form 


F  “  f  <0  (c«  Kz}  exp  (Ct  ‘  "0):l 

yields 

[3/3t  +  (V/r)  3/30]  F  =  i  M  (r)  F 
and 
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^  F  '  V  {  cos  Kz  }  •»»  0  <“  '  “» 


where 


M  (r)  =  C  -  m  V/r  (4.4) 

Alf  =  frr  +  fr/r  "  [^2/r^  +  ^  f* 

Using  these  the  following  ordinary  differential  equations  are  obtained: 


[Re-1  Urr'2)  -  i  M]  u  +  (2/r)(V+1  m  r"1  Re"1)  v  -  pp  =  0  (4.5a) 
[Re"1  (A^r'2)  -  1  H]  v  -  ( 3V/3r+V/r+2i  m  Re-1/rZ)u  +  i  m  p/r  =  0  (4.5b) 
[Re-1  Aj-i  H]  w  +  K  p  =  0  (4.5c) 
(r  u)p  -  i  m  v  +  K  r  w  5  0.  (4.5d) 


Since  V  is  obtained  from  the  finite  difference  solution  of  (2.2),  numerical 
integration  of  these  equations  is  required.  Equations  (4.5)  are  a  homo¬ 
geneous,  sixth  order  system  with  homogeneous,  two-point  boundary  conditions 
(4.2)  at  r  =  1  and  (4.3)  at  r  =  0.  The  system  is  not  self-adjoint.  Although 
there  is  no  proof,  we  assume  the  eigenvalues  form  a  denumerable,  discrete 
spectrum.  With  n  used  to  order  the  spectrum,  the  e.v.  are  C^;  for  brevity  the 

index  n  will  sometimes  be  omitted.  Eqs.  (4.5)  must  be  solved  for  Cn  and  the 
eigenfunctions  un,  vn,  wft,  pfl  given  V,  c/a,  Re,  m  and  k;  t  enters  only  through 
V. 


If  the  system  were  self-adjoint,  the  index  n  would  indicate  the  radial 
mode  number,  where  the  mode  can  be  defined  in  terms  of  the  number  of  zeros  of 
the  eigenfunctions.  For  a  system  which  is  not  self-adjoint,  it  is  unclear  how 
to  define  a  mode  in  a  general  and  unambiguous  manner.  However,  the  concept  of 
mode  is  a  very  useful  one  and  will  be  employed  here.  The  mode  can  be  identi¬ 
fied  by  calculating  the  e.v.  at  large  t,  where  the  radial  mode  number  is  known 
from  the  solid-body  rotation  solution,  and  then  tracked  as  t  decreases. 

Consider  the  inviscid  limit  of  Eqs.  (4.5),  Re"1  =  0.  The  order  of  the 
system  is  reduced  to  two,  and,  in  fact,  a  single  2nd  order  equation  is 
obtained  (Reference  3);  it  is  analogous  to  the  Rayleigh  equation  in 
hydrodynamic  stability  theory.  The  coefficient  of  the  highest  order  deriva¬ 
tive  contains  M  =  C  -  m  V/r  and  the  C  is  then  real.  If  M  *  0  has  a  root  in 
0  <  r  <  1,  the  equation  has  a  singular  point  and  the  theory  breaks  down  in  its 
neighborhood.  Let  rc  be  the  root  of  -  m  V/r.  In  analogy  with  the  Orr- 


Sommerfeld  equation  and  its  inviscid  limit,  rc  is  called  the  critical  level 

and  its  neighborhood  the  critical  layer.  A  review  of  critical  layer  dynamics 
for  shear  flows  was  given  by  Stewartson. 14  The  physical  interpretation  of 

CR  =  m  V/r 

at  r  =  rc  is  that  the  wave  frequency  is  an  integral  multiple  of  the  angular 

frequency  of  the  basic  flow,  indicating  a  resonance.  If  m  =  0  there  is  no 
critical  layer.  If  m  =  1  there  is  a  critical  layer  up  to  a  time  for  which 
rc  =  0;  after  that  a  real  rc  does  not  exist. 

For  the  viscous  case,  (4.5)  do  not  have  a  singularity  when  M  =  U. 

However,  the  critical  layer  still  plays  an  important  role  just  as  it  does  in 
the  Orr-Sommerfeld  equation.  The  r  for  which  M  =  0  is  called  a  turning 

-1/3 

point.  Lynn3  showed  that  the  thickness  of  the  critical  layer  is  0(Re“  )  in 

direct  analogy  with  the  result  for  the  Orr-Sommerfeld  equation.  The  practical 
consequence  of  the  existence  of  a  critical  layer  is  that  the  eigenfunctions 
develop  high  frequency,  large  relative  amplitude  oscillations  there.  If  the 
integration  scheme  is  not  capable  of  resolving  these,  a  solution  to  (4.5)  will 
not  be  obtained. 

There  are  two  other  properties  of  (4.5)  which  make  their  numerical 
integration  not  straightforward:  (1)  They  have  a  singularity  at  r  =  0;  (2) 

The  coefficient  of  the  highest  order  derivative.  Re'1,  is  small,  which  classi¬ 
fies  them  as  stiff  equations.  Both  of  these  properties  exist  when  the  basic 
flow  is  solid  body  rotation.  They  were  discussed  in  detail  in  Reference  1,  so 
only  a  summary  of  these  properties  and  the  methods  of  treating  than  will  be 
given. 

According  to  the  theory  of  differential  equations  (Reference  lb), 
r  =  0  is  a  singularity  of  the  second  kind.  A  singular  transformation  was 
found  that  transforms  the  system  into  one  with  a  singularity  of  the  first 
kind.  This  is  easier  to  treat,  but  the  transformed  system  has  characteristic 
values  that  are  not  all  distinct  and  differ  by  integers:  in  that  case  there 
is  no  simple,  systematic  method  of  generating  all  solutions.  (See  page  136  of 
Reference  15.)  For  the  sixth  order  system  (4.5)  there  should  be  6  independent 
solutions  near  r  =  0  and  because  of  the  boundary  conditions,  3  of  them  should 
be  singular  and  the  other  3  regular.  It  was  more  efficient  to  assume  power 
series  expansions  near  r  =  0  and  determine  the  3  regular,  linearly  Independent 
solutions.  The  only  significant  difference,  compared  to  the  solid-body  rota¬ 
tion  case  in  Reference  1,  is  that  in  the  present  case  V  must  be  expressed  as  a 


14.  K.  Stewarteon,  "Marginally  Stable  Invisoid  Flows  with  Critical  Layer'd," 
Journal  of  Applied  Mathematics,  Vol.  27,  1981,  pp.  153-176. 

15.  E.  A.  Coddington  and  N.  Levinson ,  Theory  of  Ordinary  Differential 
Equations,  McGraw-Hill  Bock  Co.,  New  York,  1956. 
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power  series  In  r.  Numerical  values  of  these  power  series  solutions  at 
r  =  e  are  used  as  initial  conditions  for  the  numerical  integration  of  (4.5) 
over  the  interval  e  <_  r  _<  1.  The  number  of  terms  in  the  series  and  e  are 
parameters  of  the  integration.  A  convergence  test  is  applied  for  each 
integration.  Formulae  and  recurrence  relations  for  the  power  series 
coefficients  are  given  in  Appendix  A. 

The  stiff  nature  of  (4.5)  was  discussed  in  Reference  1, and  ortho¬ 
normalization  was  applied  to  insure  that  the  initially  linearly  independent 
solutions  remain  so  during  the  Integration.  No  additional  discussion  is 
needed  for  the  spin-up  case.  General  rules  cannot  be  given  for  the  number  of 
subintervals,  N,  and  the  integration  step  size  which  are  required  for  a  given 
accuracy  in  the  solution  for  C;  they  depend  on  a  large  number  of  parameters, 
in  an  interactive  manner.  What  can  be  said  is  that  spin-up  and  larger  Re 
require  larger  N  and  smaller  step  size  than  solid-body  rotation  and  smaller 
Re, 


V.  SOLUTION  OF  THE  EIGENVALUE  PROBLEM 

The  procedure  for  solving  (4.5)  for  C,  and  the  eigenfunctions  u,  v,  w,  p 
is  the  same  as  that  given  in  Reference  1  for  solid-body  rotation.  The  major 
steps  in  the  procedure  are  merely  listed  here  except  for  the  "first  guess" 
strategy;  this  will  be  described  in  detail  since  it  has  no  analog  in  the  solid- 
body  rotation  case. 

For  numerical  integration,  (4,5)  are  put  into  canonical  form,  i.e.,  six 
first  order  equations  in  the  manner  of  Reference  1.  A  value  of  C  must  be 
specified,  the  first  guess.  Next,  the  power  series  solutions  are  evaluated  at 
r  =  e.  With  these  as  initial  conditions,  (4.5)  Is  integrated  for  a  specified 
N  and  integration  step  size  over  £<  r  <  1.  Orthonormalization  insures  that  3 
linearly  independent  solutions  are  generated.  These  are  tested  to  see  if  the 
boundary  conditions  are  satisfied  at  r  «  1;  this  test  requires  that  the 
characteristic  determinant,  Z  =  0;  l  Is  defined  in  Appendix  B.  If  Z  *  0  to 
within  a  certain  tolerance,  a  new  value  of  C  is  obtained  using,  essentially, 
Newton's  method.  The  iteration  then  proceeds  until  Z  *  0  and  two  successive  C 
values  differ  to  within  specified  tolerances.  The  maximum  number  of 
Iterations  allowed  is  typically  15.  For  a  convergent  Iteration  process,  the 
first  guess  must  be  reasonably  close  to  the  converged  C. 

Obtaining  a  first  guess  for  a  spin-up  e.v.  at  an  arbitrary  t  would  be 
quite  difficult.  Usually, a  spin-up  history  is  needed,  i.e.,  C  vs  t.  The 
computation  is  then  started  at  large  t,  where  a  first  guess  can  be  obtained 
from  the  C  for  solid-body  rotation  or  an  approximation  to  it.  (See  Reference  1, 
page  IB.)  For  smaller  t,  the  previously  computed  C,  or  the  two  previous  ones, 
can  be  used  to  obtain  the  first  guess.  For  some  cases,  at  early  times,  the 
required  time  interval  between  successive  calculations  of  C  becomes  small,  and 
additional  searching  for  a  first  guess  may  be  necessary  for  the  iteration 
process  to  converge. 

The  time  to  calculate  one  C  depends  mainly  on  the  first  guess  for  C,  t, 

N,  integration  step  size,  tolerances,  and  severity  of  convergence  test. 

Typical  CPU  times  for  one  C  calculation  are  one  minute  on  the  COC  7600  and  15 
minutes  on  the  VAX.  The  output  of  the  standard  program  is  only  the  e.v.,  C. 
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A  separate  program  is  used  to  print  out  the  eigenfunctions  with  an  option  to 
plot  real  and  imaginary  parts  of  w  and  p  vs  r. 

Some  results  for  eigenfunctions  and  e.v,  will  now  be  shown.  In  present¬ 
ing  these,  it  is  convenient  to  designate  the  3  wave  numbers  by  the  triplet 
(k,  n,  m)  corresponding  to  the  wave  numbers  in  the  (2,  r,  0)  directions. 


The  eigenfunctions,  of  course,  are  determined  only  to  within  a  multi¬ 
plicative  constant.  The  eigenfunctions  have  not  been  normalized  in  this 
program.  In  Figure  4  the  eigenfunction  Real  (w)  =  wR  vs  r  is  shown  at  t  = 

7000,  1000,  and  400  for  Re  =  5  x  10^,  c/a  =  2.679  (so  that  t  =  3788)  and 

mode  (6,1,1).  In  Figure  4a  t/t  =  1.85,  C„  =  8.354  x  10~2,  C,  =  8.363 
_4  s  k  1 

x  10  and  there  is  no  critical  layer*,  the  variation  of  wR  through  the 

boundary  layer  can  be  barely  discerned  on  this  scale.  In  figure  4b  t/t  = 

-2  -?  s 
.26,  Cu  =  7.060  x  10  ,  CT  =  2.269  x  10  r  =  0.44  as  indicated  by  the  arrow 

K  1  C 

in  the  figure.  The  boundary  layer  is  less  discernible  here  than  in  Figure 
4a.  In  Figure  4c  t/t  =  .11,  CR  =  3.015  x  10"2,  Cj  =  8.066  x  10'2,  and 
r  =  0.66;  away  from  the  critical  layer  wD  is  not  equal  to  zero  although  it 

C  K 

appears  to  be  on  this  scale.  For  the  conditions  of  Figure  4c,  if  n  =  628 
rad/sec  (100  Hz),  t  -  0.64  sec.  For  both  Figures  4b  and  4c,  the  critical 
layer  oscillations  are  centered  at  r  =  r  . 


Some  results  for  CR  are  presented  next  for  the  two  sample  cases  discussed 


in  the  introduction.  In  Figure  5  CR1  vs  t  is  shown  for  Case  1  with  c/a  = 
3.120,  Re  =  4  x  104 ,  mode  (3,1,1)  and  Case  2  with  c/a  =  5.200,  Re  => 


2  x  10  ,  mode  (5,1,1).  Cases  1  and 
turbulent  Ekman  layer  compatibility 


2  are  calculated  with 
conditions,  (2.8)  and 


the  laminar  and 
(2.9),  respec¬ 


tively.  For  Case  2  there  is  a  maximum  in  CR1  at  t  =  1500,  a  phenomenon  first 


determined  by  this  investigation.  If  the  calculation  were  continued  for 
smaller  times,  the  curve  would  go  through  the  origin  because  Cw  “  0  for  t  = 


0.  For  Case  1  a  maxinvuin  in  CR^  is  not  found  for  t  _>  180,  the  last  point 
computed.  The  maximum  occurs  for  t  <  130  since  CR  3  0  for  t  =  0.  For  both 


cases,  of  course,  CRj  approaches  the  result  for  solid-body  rotation  as  t  ♦  <=>. 


The  significance  of  a  maximum  in  the  CR^  vs  t  curve  is  related  to  the 

necessary  condition  for  projectile  instability  discussed  at  the  end  of  the 
introduction.  If  the  projectile  nutational  frequency  is  less  than  the  maximum 
of  CR,  there  are  two  times  at  which  instability  might  develop.  Whether  or  not 

it  does,  at  either  time,  is  not  within  the  scope  of  this  report. 

Finally,  some  results  for  another  set  of  parameters  are  presented  to  show 
the  spin-up  e.v.  histories  for  n  =  1  and  2  and  to  illustrate  the  sensitivity 
of  the  iteration  process  convergence  to  the  first  guess  of  C  that  sometimes 
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occurs  at  small  t.  These  parameters  were  chosen  by  Or.  William  P.  D'Amico  for 
a  firing  program,  proposed  a  few  years  ago,  to  be  conducted  in  the  aero¬ 
dynamics  range  for  the  purpose  of  obtaining  data  on  spin-up  instability  of 
projec.ti les.  Unfortunately,  the  program  was  not  fired, so  that  no  such  data 
exist  to  compare  with  the  theory.  However,  the  e.v.  calculations  made  for 
this  case  are  of  interest. 

In  Figure  6,  and  are  presented  as  time  histories  of  CR  and  Cj  with 

c/a  =  3.297,  Re  =  4974,  modes  (3,1,1)  and  (3,2,1).  For  this  case  the  Ekman 
layers  would  be  laminar,  so  that  (2.8)  is  used  in  the  spin-up  basic  flow 
calculation.  The  curve  has  an  unusual  feature.  For  180  <  t  <  220  it  has 

a  shallow  minimum  and  maximum;  on  the  scale  of  this  figure  appears  to  be 

constant  over  this  range.  The  variation  of  Cj^  is  rather  typical.  Obtaining 

a  first  guess  for  C,  for  which  the  iteration  process  would  converge,  was  made 
more  difficult  because  was  not  monotonic.  However,  the  calculation 

of  was,  relative  to  C^,  considerably  more  difficult  starting  at  about  t  = 

250.  The  range  of  values  of  CR^  »  C j ^  which  a  first  guess  for  had  to 

lie,  in  order  for  the  iteration  process  to  converge  to  the  n  =  2  root  of  Z  s 
0,  steadily  decreased.  This  is  shown  in  the  Cj  vs  CR  plot  in  Figure  7.  The 

curve  is  the  e.v.  trajectory  for  180  <  t  <  250;  the  areas  in  which  a  first 
guess  must  lie  are  shown  for  four  values  of  t.  As  indicated  for  t  =  200  and 
250,  the  areas  are  off-scale.  For  t  a  180  the  rectangle  has  dimensions 
(.0008,  .0007).  The  search  for  a  first  guess  must  be  done  with  differences 
less  than  .001.  The  difficulty  of  obtaining  a  first  guess  is  compounded  by 
the  fact  that  CI2  is  nonmonotonic  near  t  a  200.  When  the  iteration  process 

did  not  converge  to  the  n  =  2  root  of  2  =  0,  it  almost  always  converged  to  the 
n  =  1  root;  sometimes  It  converged  to  a  root  for  larger  n,  and  a  few  times  It 
did  not  converge.  There  was  a  very  predominant  preference  for  the  n  »  1  root. 

To  visualize  Z  (C),  |Z|  (CR,  Cj)  is  plotted  in  Figure  8.  Only  the 

shaded  portion  is  the  jZ|  surface;  the  planar  regions  from  the  surface  to  |Z| 

=  0  are  artifacts  of  the  plotting  program.  The  surface  is  plotted  in  the 
neighborhood  of  the  n  »  2  root  for  t  °  180,  viz.,  CR^  B  .101673,  » 

.034867.  Only  some  of  the  values  on  the  coordinate  lines  are  labeled.  The 
root  lies  at  the  bottom  of  the  steep  valley  (not  visible  in  the  plot).  To  the 
right  of  it  is  a  steep  peak,  and  this  helps  to  explain  why  most  guesses  did  not 
converge  to  the  n  »  2  result.  In  a  sense  the  valley  is  shaded  by  the  peak. 
This  case  illustrates  the  fact  that  for  early  times  the  solution  to  the  e.v. 
problem  may  become  more  laborious  than  for  larger  times. 


VI.  DISCUSSION 

In  this  report  the  theory  and  a  method  for  the  solution  of  the  spin-up 
e.v.  problem  are  given.  The  theory  for  the  perturbation  flow  is  a  linear  one 
but  not  for  the  basic  flow.  Because  viscous  effects  are  important  In  the  e.v. 
problem  at  the  cylinder  wall  and  critical  layer,  viscous  perturbations  are 
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used.  Several  assumptions  are  made,  and  some  limitations  are  contained  in  the 
theory.  The  theory  and  methods  are  successful  in  the  sense  that  they  provide 
results,  as  shown  in  the  previous  section,  that  are  physically  meaningful  and 
do  not  violate  intuition  or  the  "physics  of  the  problem."  Other  investigators 
have  worked  on  this  problem  without  success  in  that  sense.  The  more  important 
gauge  of  success  is  validation  by  comparison  with  either  numerical  simulation 
or  experimental  results. 

For  the  solid-body  rotation  case  and  m  =  0,  a  detailed  validation  was 
given  in  Reference  7  using  finite  difference  solutions  to  Navier-Stokes 
equations.  This  provides  confidence  in  the  treatment  of  (4.5),  and  the  solid- 
body  rotation  endwall  correction,  for  that  case.  Spin-up  and  m  =  0  could  be 
validated  in  the  same  way;  it  would  require,  relatively,  more  analysis  to 
reduce  the  numerical  data.  At  the  present  time  there  is  no  numerical 
simulation  to  validate  the  m  -  1  case. 

Although  the  importance  of  the  spin-up  problem  for  liquid-filled  pro¬ 
jectiles  has  been  recognized  since  the  early  1960's,  very  few  experiments 
relevant  to  spin-up  have  been  conducted  at  BRL,  and  none  of  these  were 
directed  to  the  modes  of  oscillation.  The  only  published  results  seem  to  be 
those  quoted  in  References  4  and  10.  Experiments  in  a  spin  generator 
determined  particle  paths  during  spin-up,  and  range  tests  provided  spin  decay 
histories;  these  were  used  by  Wedemeyer10  to  validate  his  theory  for  the  spin- 
up  basic  flow,  described  in  section  II.  Further  validation  of  that  flow  was 
obtained  with  measurements  of  velocity  using  the  LDV,  e.g..  Reference  16. 

The  only  attempts  to  measure  the  e.v.  were  those  of  Aldridge17  and  the 
further  development  of  that  work  reported  in  Reference  6.  In  Reference  17, 
measurements  of  CR  (not  Cj)  are  given  for  modes  (2,1,0),  (2,2,0),  and  (4,2,0) 

over  the  range  1.4  <  t/ t^  <  ®  and  compared  with  Lynn's  calculation  for 

inviscid  perturbations  and  m  =  0.  For  larger  times  the  difference  was  small; 
but  the  trend,  as  t  decreases,  was  for  the  differences  to  grow.  Later, 
Aldridge  found  that  a  systematic  error  existed  in  the  analysis  of  the  data. 
Some  of  these  experiments  were  redone  using  a  different  data  analysis 
technique  and  reported  in  Reference  6;  they  are  presented  below.  For  in  =  0 
there  is  no  critical  layer  and  t/t$  *  1.4  is  larger  than  one  would  want  for  a 

spin-up  test.  In  Reference  6  measurements  of  and  Cj  are  given  for  modes 

(2,1,0),  (1,2,1)  and  (1,1,1).  The  first  of  these  replaces  those  of  Reference 
17.  The  latter  two  were  excited  using  a  rotating  cylinder  with  a  precessing 
top.  In  Table  1,  these  experimental  results  are  presented  together  with 
calculated  results  using  the  method  described  here.  The  smallest  t/t$  = 


1 (i.  U.  8.  k'atkina  and  G.  8.  Hueaey ,  "Spin-Up  ft*om  Root  in  a  Cylinder,"  The 
PhyaieB  of  Fluids,  Vol.  SO,  Ho,  10,  Part  1,  1977,  pp,  1SS6-1604 . 

17.  K.  D.  Aldridge,  "Experimental  Verification  of  the  Inertial  Oaoillatione 
of  a  Fluid  in  a  Cylinder  During  Spin-Up BRL  Contmot  Report  Ho  $73, 
Alwrdeen  Proving  Ground,  Maryland,  September  197S  (AD  AO  18797).  Alao 
Geophya.  Aotrophye.  Fluid  Dynamics,  Vol .  R,  1977,  pp.  279-301. 
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1.51.  The  experimental  CR  has  an  estimated  error  of  less  than  1%;  the  error 

In  Cj  can  be  as  much  as  10%.  The  differences  between  the  experimental  and 

calculated  frequencies  are  less  than  2%.  For  the  damping  the  difference  is  a 
factor  of  3  for  t/t$  «  1.51;  this  could  be  expected  because  we  have  neglected 

the  Ekman  layers,  as  discussed  in  Sections  III  and  IV.  Solid-body  rotation  is 
achieved,  essentially,  for  t/t$  *  4.75;  thus,  the  method  of  Reference  1  can  be 

used,  including  the  solid-body  rotation  endwall  correction  {not  Ekman  layer 
correction).  This  gives  the  Cj  in  the  last  column  of  Table  1,  and  the 
difference  is  now  about  10%. 


TABLE  1.  EXPERIMENTAL  AND  CALCULATED  EIGENVALUES 
(2,1,0)  MOOE  •  SPIN-UP 
Re  ~  43,000  c/a  =  0.995 


Frequency 

Damping 

t/ts 

Exp. 

Calc. 

Exp. 

Calc.  Uncorr. 

Calc.  Corr 

4.75 

1.269 

1.268 

.0061 

.00297 

.007 

3.04 

1.266 

1.265 

.0062 

.00295 

2.02 

1.236 

1.242 

.0068 

.00311 

1.68 

1.202 

1.219 

.0101 

.00334 

1.51 

1.184 

1.202 

.0105 

.00342 

The  comparisons  of  CR  for  the  ra  =>  1  modes  are  given  in  Table  2.  For  each 

t.  Re  varies  slightly;  the  values  given  are  representative.  Again,  the 
differences  between  the  experimental  CR  and  our  calculated  results  are,  at  most, 

2%.  Further  validation  is  needed  for  smaller  t/ t s .  for  both  smaller  and 

larger  Re  and  for  larger  c/a. 

For  projectile  applications  m  *  1  is  the  pertinent  azimuthal  wave  number 
and,  to  apply  the  necessary  condition  for  projectile  instability,  CR  is 

required.  The  comparison  shown  in  Table  2  gives  the  significant  result  that 
the  method  presented  here  is  validated  within  the  limitations  noted. 

The  lack  of  projectile  tests  in  ballistic  ranges  has  already  been  noted. 
Field  tests  wherein  data  is  obtained  from  yawsondes  can  provide  sane  informa¬ 
tion,  but  not  as  detailed  as  range  tests.  For  the  field  tests  on  155mm 
projectiles  mentioned  in  the  introduction,  some  of  the  parameters  were  chosen 
on  the  basis  of  our  calculations  to  obtain  information  on  instabilities  during 
spin-up.  The  results  are  in  publication  by  D'Amico.16 


IB. 


V.  P.  D'Amico,  " Flight  Data  on  Liquid-Pit  1 sd  Shell  for  Spin-Dp 
Instabilities,"  BRL  Report  in  prepaixition,  193$. 


TABLE  2,  EXPERIMENTAL  AND  CALCULATED  E IGENFREQUENC I E S 

c/a  =  0.60C 


Mode  (1,2,1) 

t/ts 

Re  s  50,000 

Calculated 

U-Cr) 

Experimental 

(i-cR) 

Di fference 

1,473 

.8474 

.8658 

2.1% 

1.804 

.8352 

.8488 

1.6% 

2.800 

.8276 

.8305 

.35% 

4.947 

.8255 

.8295 

.49% 

Mode  (1,1,1) 

Re  ~  30,000 

1.357 

1.3359 

1.3487 

.95% 

1,634 

1.3707 

1.3736 

.2% 

3.780 

1.4275 

1.426 

.1% 

4. HO 

1.4278 

1.4290 

.1% 

Free  oscillation  of  the  fluid  during  spin-up  was  considered  in  this 
report.  For  projectile  motion  the  forced  oscillation  problem  must  be 
sol  .’ed.  The  modal  solutions  to  (4.5)  would  be  needed  there  also,  but  the 
boundary  conditions  are  nonhomogeneous.  As  shown  in  Reference  19  for  the 
solid-body  rotation  case,  the  same  techniques  discussed  in  this  report  are 
needed  in  the  forced  oscillations. 
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V  vs  r  for  Case  1  with  c/a  3  3.120,  Re  =  4  x  10  at 
Three  Values  of  t;  t  *  1245.  Laminar  Ekman  Layer. 


28 


I 

.0 


figure  2 


i  i  I  i  I 

.2  ,4  .6  .8  1.0 


r 


V  vs  r  for  Case  2  with  c/a  =  5.200,  Re  K  2  x  10^  at 
Three  Values  of  t;  a  2700.  Turbulent  Ekman  layer. 
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V*  vs  z  +  c/a  =  z'  for  c/a  =  1,  Re  =  10  ,  tg  a  <-0° 

at  (a)  t  62.5,  (b)  t  =  250.  The  Solution  is  Obtained 
From  the  Program  in  Reference  11. 
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Figure  7.  Trajectory  of  C2  With  t  as  Parameter;  c/a  =  3.297,  Re  =  4,974,  Mode  (3,2,1). 

Approximate  Regions  in  Which  First  Guess  for  C2  Must  Lie  to  Have  Convergence 
to  n  =  2  Mode  Are  Shown  for  Four  Times. 


Surface  |Z|  (CR,  Cj)  for  c/a  =  3.297  and 
he  Neighborhood  of  the  Root  CR?  =  .101673 
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APPENDIX  A.  POWER  SERIES  COEFFICIENTS 


The  need  for  the  power  series  expansions  near  r  =  0  of  u,  v,  w,  and  p  is 
described  in  Section  IV.  The  coefficients  of  these  series  are  recorded  here; 
they  are  obtained  by  substituting  the  series  Into  (4.5).  V  must  also  be 
expressed  in  a  power  series 

V  ■  E  V.  r> 
j  -  1  J 

which  is  permissible  since  V  is  an  analytic  function  of  r,  near  r  =  0,  except 
that  it  may  not  be  uniformly  so  for  t  -*•  0.  The  coefficients  Vi  could  be 

u 

determined  from  (2.2)  in  terms  of  aV/3r,  3  (3V/3t)/3r,  etc.,  at  r  =  0.  These 
would  have  to  be  found  numerically  from  the  solution  V  (r,t).  Instead  we 
evaluate  the  V.-  directly,  fitting  a  quartic  to  V  (r,t). 

Because  the  boundary  conditions  (4.3)  are  different  for  m  ®  0  and  m  =  1, 
the  coefficients  must  be  given  separately  for  these  two  cases. 

m  =  0: 

The  forms  of  the  power  series  are 

u  =  L  u.  rJ  v  »  £  v,  r* 

j  *  1  J  j  »  1  J 

"  ■  "o  4  .£  “i  rJ  P  "  P0  *  .1,  i>j  ^ 

j  a  1  J  J  »  1  J 

The  coefficients  are  given  by  the  following  sequence  of  formulas,  where 
wq,  pQ,  and  v^  are  arbitrary  complex  constants: 

H  =  X2  Re'1  +  i  C 

u1  =  -  (K/2)wq,  w,  «  0 

vl  a  V  Pi  *  0 

u2  «  0,  w?  «  (Re/4)  (Kw0  -  Kp0) 

=*0,  p2  •  (1/2)  Ku1  ♦  (K/Rc)  w2  -  Vj  v1 

For  j  *  1,  2,  3  ... 

wj+2  “  ke  U+2>’2 
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uj+2  =  “K^+3^  1  wj+i 

2  1  J  ♦  1 

vJ>2  =  Re  [(j*  +  4j  ♦  3)]'1  [Hvj  ♦  £  (Ul)  uj+1-£] 

i  **  1 

Pjt2  ■  (j+2)"‘  W*)'1  (H.j  -  Kpj)  -  HUj+1  +  2  ^  V,  vj+2_£] 

m  =  1: 

The  solution  for  the  azimuthal  mode  m  =  1  is  given  by 


■i?o  »/• 

*  ■  j?0  bjr'1 

■  i=oVj' 

p  =  .  E  e  .  r3 

J  =  0  j 

Summations  are  taken  from  j  =  0  for  convenience. 

The  coefficients  are  given  by  the  following  sequence  of  formulas,  where 
bQ>  dp  and  e^  are  arbitrary  complex  constants: 


=  ib  , 

0* 

b  a  b  , 

0  o 

do  .  0, 

eo=° 

-  0, 

bl  -  0, 

di  "  V 

e^  a  e 

H  =  Rq'1  K2  +  i  C 

a2  -  -(l/4)Kd1  ♦  (1/8)  Re  [bQ  (iH  -  V  )  ♦ 

b2  -  -i( l/4>Kd1  -  (31/8)  Re  (bQ  (iH  -  Vj)  +  e^ 
For  j  «  0,  1,  2,  ... 


Re  [Hajn  - 
Re  CHbj+1  + 


j  *  1 

E 


8  =  0 


Vt.l  *  ‘'j,!.,)! 

Vi 
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Nj  =  i\.  -  i{j+2)N2j]/[(j+l)  (j+3)J 

dj+2  -  Re  [Hd.  -  tej  -  i  45o  V£+1dj.,]/C(j+l)  ( J+3)] 

aj+3  =  -[N3  +  (j+4)  Kd.+2]/(j2  +  8j  +  15) 

bj+3  =  i[Kdj+2  +  (j+4)Nj]/(j2  +  8j  +  15) 

ej+2  =  C2aj+3  +  1  (J2  +  6j  +  7)  bj+3  -  iN2j]/Re 
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APPENDIX  B.  CHARACTERISTIC  DETERMINANT,  Z 

The  characteristic  determinant,  Z  (C),  denoted  by  D  (C)  In  (25)  of 
Reference  1,  is  defined  here  for  convenience.  In  the  integration  procedure  we 
calculate  three  linearly  Independent  solutions  that  satisfy  three  boundary 

conditions  at  r  =  0.  These  solutions  are  denoted  by  (u.,  u.  -  i  v.,  w.), 

J  J  J  J 

j  =  1,  2,  3;  pj  is  not  needed  here.  The  no-slip  conditions  at  the  sidewall 

require  the  velocity  components  to  vanish  at  r  =  1.  These  components  are 
linear  combinations  of  the  three  solutions  above.  Thus 

Y1  1^(1)  +  y2  [u1(1)  -  i  Vjd)]  +  y3  Wjd)  =  0 

Y1  +  y2  !-u2^  "  1  +  y3  W2  W  =  0  (B1) 

Yj  u3(l)  +  y2  [u3(1)  -  i  v3( 1)]  +  y3  w3(1)  =  0, 

where  y^  y2»  and  y3  are  constants.  In  order  that  a  non-trivlal  solution  for 
Yji  y2.  and  y3  be  found,  the  determinant  of  the  coefficients  in  (Bl),  which  is 
designated  Z,  must  equal  zero. 


LIST  OF  SYMBOLS 


n 

N 


P 

P* 

P 

P 


o 


rl 

P* 


r 


Re 

t 


cross-sectional  radius  of  cylinder 

half-height  of  cylinder 

a  Cg  +  i  Cj,  complex  eigenvalue 

disturbance  decay  rate/« 

disturbance  frequency/ ft 

C^'s  for  n  =  1  and  n  =  2  modes,  respectively 

=  v/ilc2,  Ekman  number 

axial  mode  number  (see  (4.1)) 

=  k  (a/c)  Re”1/f2  (see  (2.8)) 

=  0.035  (a/c)  Re~1/5  (see  (2.9)) 

=  k  ita/(2c)  (see  (4.1)) 
azi.wuthal  mode  number  (see  (4.1)) 

=  C  -  mV/r 
radial  mode  number 

number  of  subintervals  used  in  orthonorraalization 
function  describing  radial  variation  of  p'  (see  (4.1)) 

3-D  pressure  perturbation/(pa2a2) 

pressure/(pa2n2)  of  approximate  basic  flow  (see  (2.3)) 

pressure/ ( pa 2q2)  in  3-D  Navier-Stokes  equations  (see  (3.1)) 

zeroth  order  approximation,  in  l/tg,  to  P*  (see  (2,6)) 

coefficient  in  first  order  approximation  to  P*  in  (2,6) 

pressure/(pa2Q2)  of  exact  basic  flow  (see  (2.1)) 

radial  coordinate  /  a 

critical  level  radius  /  a 

=  a 2a/v.  Reynolds  number 

time  x  Q 

1  111  ■  18  — m  ii  , 
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characteristic  flight  time  of  projectile 

t  =  (2c/a)  Re1//2,  t$  =  t g/ft  =  spin-up  time  (laminar) 

t$t  =  (28.6c/a)  Re1//5,  =  t$t/8  =  sP''n,-uP  time  (turbulent) 

t*  2t/t$ 

u,  v,  w  functions  describing  radial  variation  of  u' ,  v‘,  w'  (see  (4.1)) 

u1,  v‘,  w*  radial,  azimuthal,  axial  velocity  components  x  1/ ( a  n)  of  3-D 

perturbed  flow  (see  (3.2)) 

U,  V,  W  radial,  azimuthal,  axial  velocity  components  x  1/ ( a  ft)  of 

approximate  basic  flow  (see  (2.2)) 

U . ,  V,,  W.  coefficients  in  first  order  approximations,  in  i/t.  to  U*,  V*, 

1  1  1  W*  (see  (2.6))  s 

Uw,  V  radial  and  azimuthal  velocity  components  x  1/ ( a  ft)  of  basic  flow 

*  without  diffusion  given  by  (2.4) 

U,  V,  W  radial,  azimuthal,  axial  velocity  components  x  1/ ( a  ft)  of  the 

3-D  Navier-Stokes  equations  (see  (3.1)) 

U*,  V*,  W*  radial,  azimuthal,  axial  velocity  components  x  1/ ( a  ft)  of 
exact  basic  flow  (see  (2.1)) 

VQ  zeroth  order  approximation,  in  l/ts,  to  V*  (see  (2.6)) 

z  axial  coordinate/a  (z  3  0  at  cylinder  midplane) 

z*  «  z  +  c/a 

Z  characteristic  determinant  whose  roots  are  eigenvalues 

e  value  of  r  where  numerical  integration  is  initiated 

o  azimuthal  angle 

x  constant  in  expression  for  radial  velocity  in  approximate  basic 

flow  with  laminar  Ekman  layer  (see  (2.8)) 

v  kinematic  viscosity  of  fluid 

p  density  of  fluid 

ft  spin  rate  of  cylinder 
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